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Abstract 


We investigate the mass spectra of closed-charm mesons with = 
1 , for hybrid charmonium, molecules, and diquark-antidiquark op¬ 

erators, in quenched lattice QCD with exact chiral symmetry. For 
two lattice volumes 24^ X 48 and 20^ X 40, each of 100 gauge con¬ 
figurations generated with single-plaquette action at /3 = 6.1, we com¬ 
pute point-to-point quark propagators and measure the time-correlation 
functions of these exotic meson operators. For the molecular operator 
{(q757jc)(c75q) — (c757jq)(q75c)}, it detects a resonance with mass 
around 4238±31 MeV, which is naturally identified with y(4260). Fur¬ 
ther, for any molecular and diquark-antidiquark operator, it detects 
heavier exotic charmed mesons, with (cscs) around 4450 ± 100 MeV, 
and (cccc) around 6400 ± 50 MeV. 


PACS numbers: 11.15.Ha, 11.30.Rd, 12.38.Gc 
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1 Introduction 


Recently, a new state y(4260) with a mass of 4259(8)(4) MeV and a width 
88(23)(5) MeV, in association with an initial state radiation photon, has been 
observed by BaBar collaboration in annihilation [1], This immediately 
implies that its = 1 . From the experimental and theoretical spectrum 

of (cc) states, V(4260) can hardly be interpreted as one of the radial and/or 
orbital excitations of (cc). Thus it is most likely an exotic (non-gg) meson. 

So far, theoretical models to understand Y (4260) are: (i) hybrid charmo- 
nium ccg [2, 3], (ii) P-wave excitation of the diquark-antidiquark [cs][cs] [4]; 
and (iii) molecule composed of 2 mesons, e.g., pxci [5]. 

Now the important question is whether there is a resonance around 4260 
MeV with = 1 , in the spectrum of QCD. In this paper, we investigate 

the mass spectra of several interpolating operators whose lowest-lying states 
have = 1 , in lattice QCD with exact chiral symmetry [6, 7, 8, 9, 10]. 

For two lattice volumes 24^ x 48 and 20^ x 40, each of 100 gauge conhgurations 
generated with single-plaquette action at /? = 6.1, we compute point-to-point 
quark propagators for 30 quark masses in the range 0.03 < m^a < 0.80, and 
measure the time-correlation functions of the exotic meson operators which 
can overlap with V(4260). The inverse lattice spacing a~^ is determined with 
the experimental input of /tt, while the strange quark bare mass m^a = 0.08, 
and the charm quark bare mass rrica = 0.80 are hxed such that the masses 
of the corresponding vector mesons are in good agreement with 0(1020) and 
7/0(3097) respectively [11]. Our scheme of computing quark propagators has 
been outlined in [12]. 

Note that we are working in the quenched approximation which in principle 
is unphysical. However, our previous results on charmed baryon masses, and 
also charmed meson masses and decay constants (theoretical predictions) [11] 
turn out to be in good agreement with the experimental values. This seems to 
suggest that it is plausible to use the quenched lattice QCD with exact chiral 
symmetry to investigate the mass spectra of the charmed meson operators con¬ 
structed in this paper, as a hrst step toward the unquenched calculations. The 
systematic error due to quenching can only be determined after we can repeat 
the same calculation with unquenched gauge conhgurations. However, the 
Monte Carlo simulation of unquenched gauge conhgurations for lattice QCD 
with exact chiral symmetry, on the lattices 20^ x 40 and 24^ x 48 at 0 = 6.1, 
still remains a great challenge to the lattice community. Thus, in this paper, 
we proceed with the quenched approximation, assuming that the quenching 
error does not change our conclusions dramatically, in view of the good agree¬ 
ment between our previous quenched mass spectra of charmed hadrons [11] 
and the experimental values. 
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2 The Hybrid Charmonium ccg 

The local interpolating operators for hybrid charmonium with = 1 
can be constructed as cF^iC, eijkC'j^FjkC, and eijkC'j^'j^FjkC. Here the matrix¬ 
valued gluon held tensor F^y{x) can be obtained from the four plaquettes 
surrounding x on the (/i, z>) plane, i.e., 

ga^F^^{x) ~ ^[P^^{x) + P^u{x - A) + Pkiu{,x - P) + P^y{x - jl - z>) 

Ol 

- PlA^ -h) - PlA^ -^)-PlAx-p- ^)] 


where 


Pfiuix) = U^{x)Uu{x + fi)Ul{x + P)Ul{x) . 
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Figure 1: (a) The time-correlation function C{t) of the hybrid meson operator 
^ijk^'JbPjkC with = 1 , on the 24^ x 48 lattice at (3 = 6.1. The solid 

line is the double hyperbolic-cosine ht for t G [1,24]. (b) The effective mass 
Meff{t) = ln[C'(t)/C'(t-|- 1)] of C{t) in Fig. la. 


In Fig. 1, the time-correlation function C{t) of the hybrid charmonium 
operator eijkC'y^FjkC is plotted, together with the effective mass ln[C(t)/C'(t-|- 
1)], for 100 gauge conhgurations generated with single-plaquette action on 
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24^ X 48 at /3 = 6.1. Here C{t) has been averaged over z = 1, 2, 3, where in each 
case, the “forward-propagator” Ci{t) and “backward-propagator” Ci{T — t) 
are averaged to increase the statistics. The same strategy is applied to all 
time-correlation functions in this paper. The solid line in (a) is the double 
hyperbolic-cosine £t 

for t G [1,24]. It gives mi = 2977(28) MeV and m 2 = 4501(178) MeV with 
X^/d.o.f. = 0.95. In this paper, we have adopted the procedure outlined in 
Ref. [13] to perform our data htting and error estimation. Now we identity 
the lowest-lying state with since the hybrid charmonium operator with 

1 also overlaps with J/4’. Then the hrst excited state with mass 4501(178) 
MeV is identihed with the lowest-lying hybrid charmonium with 1 . Note 

that the error of the mass of the first excited state is relatively large, since it 
is obtained by double hyperbolic-cosine fit. Evidently the mass of the lowest- 
lying hybrid charmonium state with 1 is higher than 4260 MeV. Thus it is 
unlikely to be identihed with Y(4260), even though we could not rule out such 
a possibility, due to the large error bar. Nevertheless, we will report a more 
precise determination of the mass of the hybrid charmonium with 1 in a 
future publication. 


3 The Molecular Operators 


In this section, we construct three molecular operators with quark content 
(cqcq) such that the lowest-lying state of each operator has = 1 . Then 

we compute the time correlation function of each operator, and extract the 
mass of its lowest-lying state. Explicitly, these molcular operators are: 


Ml 

M 2 

M 3 


1 


{(q 7 ic)(cq) + (c 7 iq)(qc)} 


V2 

(c7ic)(qq) 

^ {(q757ic)(c75q) - (c757iq)(q75c)} 


( 1 ) 

( 2 ) 

( 3 ) 


The time-correlation function^ 


CM{t) = 

X 

is measured for each gauge conhguration, and its average over all gauge con- 
hgurations is htted to the usual formula 

_fg—mat _|_ ma(T— 

_2ma 

^Here we have neglected the cc and qq annihilation diagrams such that C(t) does not 
overlap with any conventional meson (cc or qq) states. 
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Figure 2: The ratio of spectral weights of the lowest-lying state of the molecular 
operator Mi, for 20^ x40 and 24^ x48 lattices at /3 = 6.1. The upper-horizontal 
line R = (24/20)^ = 1.728, is the signature of 2-particle scattering state, while 
the lower-horizontal line i? = 1.0 is the signature of a resonance. 


to extract the mass ma of the lowest-lying state and its spectral weight 



Theoretically, if this state is a genuine resonance, then its mass ma and spectral 
weight W should be almost constant for any lattices with the same lattice 
spacing. On the other hand, if it is a 2-particle scattering state, then its 
mass ma is sensitive to the lattice volume, and its spectral weight is inversely 
proportional to the spatial volume for lattices with the same lattice spacing. 
In the following, we shall use the ratio of the spectral weights on two spatial 
volumes 20^ and 24^ with the same lattice spacing (/3 = 6.1) to discriminate 
whether any hadronic state under investigation is a resonance or not. 

In Fig. 2, the ratio (i? = W 20 /W 24 ) of spectral weights of the lowest-lying 
state extracted from the time-correlation function of Mi on the 20^ x 40 and 
24^ X 48 lattices is plotted versus the quark mass rUgO G [0.03, 0.8]. Evidently, 
R ~ 1.0 for mqa > 0.05, which implies that there exist 1 resonances with 
quark contents (cccc) and (cscs). On the other hand, as > m„, R begins 
to deviate from 1.0 with large error. We suspect that this might be due to 
quenching artifacts as —*• m„, mostly from the scalar meson component 
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Figure 3: The mass of the lowest-lying state of Mi versus the quark mass 
niqa, on the 24^ x 48 lattice at /3 = 6.1. The solid line is the linear ht with 
xVd.o.f. = 0.34. 


(cq) or (qc) in the molecular operator Mi. Thus R could be consistent with 
1.0 if one incorporates internal quark loops, with higher statistics and larger 
volumes. If this is the case, a resonance might also exist for q = u. In the 
following, we assume this is the case, and chirally extrapolate the mass of the 
molecule Mi to the limit rriq —> m„. 

In Fig. 3, the mass of the lowest-lying state extracted from the molecular 
operator Mi is plotted versus m^a, which can be htted by the linear function 
m = Co -|- Ciuiq with /d.o.f. = 0.34. In the limit rriq —*• rUu, it gives m = 
4350(69) MeV which is compatible with the mass of F(4260). 

For rriq = rris, and rriq = rric, the time-correlation functions and effective 
masses of Mi are plotted in Fig. 4 and Fig. 5 respectively. The masses of 
the lowest-lying states are: m[(s 7 jc)(cs) -|- (c 7 js)(sc)] = 4546(30) MeV with 
X^/d.o.f. = 0.97, and m[(c 7 jc)(cc)] = 6411(25) MeV with /d.o.f. = 0.61. 

Next we turn to the molecular operator M 2 . Obviously, it must suffer 
severely from the quenched artifacts as rriq —>• rriu, due to the presence of the 
scalar meson operator (qq). Thus, we observe that in the limit rriq —>■ rriu, 
the time-correlation function becomes very noisy, and we could not obtain any 
reliable results of its masses and spectral weights. However, at rriq = rris = 
0.08a“^, we can still extract its mass and spectral weight. For the 24^ x 48 
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Figure 4: (a) The time-correlation function C{t) of the lowest-lying state of 
Ml for rriq = rus = 0.08a“^, on the 24^ x 48 lattice at /3 = 6.1. The solid line is 
the hyperbolic-cosine £t for t G [9,13] with x^/d.o.f. = 0.97. (b) The effective 
mass Meff{t) = \n[C{t)/C{t + 1)] of C{t) in Fig. 4a. 


lattice (Fig. 6), the mass and spectral weight are: m = 4581(96) MeV, and 
W = 0.71(41) X 10“^; while for the 20^ x 40 lattice, m = 4637(156) MeV, and 
W = 0.92(80) X 10“’^. Thus the molecular operator (c 7 jc)(ss) seems to detect 
a resonance with mass 4581(96) MeV (Fig. 6), which is compatible to that 
obtained from Mi with ruq = (Fig. 4). 

Next we consider the molecular operator M 3 . It is expected to suffer 
quenched artifacts less than those of Mi and M 2 , since it is composed of 
a pseudoscalar meson operator (cysg) times a pseudovector meson operator 
(q 757 jc), without any scalar meson operators like (qq) or (cq). In Fig. 7, the 
ratio (i? = W 2 Q/W 2 i) of spectral weights of the lowest-lying state extracted 
from the time-correlation function of M 3 on the 20^ x 40 and 24^ x 48 lattices 
is plotted versus the quark mass m^a G [0.03,0.8]. Evidently, R ~ 1.0 for the 
entire range of quark masses, which implies that there exist resonances with 
quark contents (cccc), (cscs), and (cucu), even though R slightly deviates 
from one at small quark masses. 

In Fig. 8 , the mass of the lowest-lying state extracted from the molecular 
operator M 3 is plotted versus m^a, which can be htted by the linear function 
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Figure 5: (a) The time-correlation function C{t) of the lowest-lying state of 
the molecular operator (c 7 jc)(cc), on the 24^ x 48 lattice at /3 = 6.1. The solid 
line is the hyperbolic-cosine £t for t G [11,17] with /d-o.f. = 0.61. (b) The 
effective mass Meff{t) = \n[C{t)/C{t+ 1)] of C{t) in Fig. 5a. 


m = Co -|- CiTUq with x^/d-o.f. = 0.25. In the limit niq —> it gives m = 

4238(31) MeV, which is in good agreement with the mass of Y (4260). 

For rriq = rUg = 0.08a“^, the time-correlation function and effective mass 
of M 3 are plotted in Fig. 9. The mass of the lowest-lying state is 4405(31) 
MeV, which is compatible to those obtained from Mi and M 2 with ruq = m^. 

4 The Diquark-Antidiquark Operator 

We construct the diquark-antidiquark operator with = 1 as 

Ya{x) = {[(f Cx5lic\xa[^CxbC^]xa + [qCysyfc^]a,a[q^C'yscj^a} (4) 

where C is the charge conjugation operator satisfying = — 7 J and 

(^^ 75 )^ = —Cxb- Here the diquark operator [q^FQIj-a for any Dirac matrix F 
is dehned as 




/3qa;/3c) 


( 5 ) 
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(a) (b) 

Figure 6 : (a) The time-correlation function C{t) of the lowest-lying state of 
(c 7 jc)(ss), on the 24^ x 48 lattice at (3 = 6.1. The solid line is the hyperbolic- 
cosine ht for t G [9,14] with x^/d-o.f. = 0.56. (b) The effective mass Meff(t) = 
\n[C{t)/C{t + 1)] of C{t) in Fig. 6 a. 


where x, {a,b,c} and {a,/?} denote the lattice site, color, and Dirac indices 
respectively, and eabc is the completely antisymmetric tensor. Thus the diquark 
(5) transforms like color anti-triplet. For F = Cys, it transforms like = O’*", 
while for F = C'y^'yi {i = 1, 2, 3), it transforms like 1“. In the limit q = Q, the 
diquark operator is replaced by (q^Fq)^;^ = eatcqxab^apqxpc- 

In Fig. 10, the ratio {R = W 2 Q/W 24 ) of spectral weights of the lowest-lying 
state extracted from the time-correlation function of 14 on the 20^ x 40 and 
24^ X 48 lattices is plotted versus the quark mass m^a G [0.03, 0.8]. Evidently, 
R ~ 1.0 for niqa > 0.05, in particular, for rriq = ms = 0.08a“^, and = 
rUc = 0.8a“^. This implies that that there exist 1 resonances with quark 
contents (cccc) and (cscs). On the other hand, as > m„, R begins to 
deviate from 1.0 with large error. This seems to suggest that the diquark- 
antidiquark operator F 4 has little overlap with the resonance detected by the 
molecular operator M 3 as mq —> However, it is unclear whether i? ~ 1 

if one incorporates internal quark loops, and with larger volumes and higher 
statistics. In the following, we assume this is the case, and obtain its mass by 
chiral extrapolation. 
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Figure 7: The ratio of spectral weights of the lowest-lying state of the molecular 
operator M 3 , for 20^ x 40 and 24^ x 48 lattices at /3 = 6.1. The upper-horizontal 
line R = (24/20)^ = 1.728, is the signature of 2-particle scattering state, while 
the lower-horizontal line i? = 1.0 is the signature of a resonance. 
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In Fig. 11, the mass of the lowest-lying state of the diquark-antidiquark 
operator I 4 is plotted versus m^a, which can be fitted by the linear function 
m = Co -|- Ciuiq with x^/d-o.f. = 0.49. In the limit rriq niu, it gives m = 
4267(68) MeV, which is in good agreement with the mass of T(4260). 

For niq = rris = 0.08a“^, and rriq = rric = 0.80a“^, the time-correlation 
functions and effective masses of the diquark-antidiquark operator are plotted 
in Fig. 12, and Fig. 13 respectively. The masses of the lowest-lying states 
are: m{[s^C 757 jc] [sCysC^] -|- [ 8 ( 7757 ^c^] [s^C 75 c]} = 4449(40) MeV with 
X^/d.o.f. = 1.03, and m{(c^C 757 ic)(cC 75 C^) + ( 0 ^ 757 ^ 0 '^)( 0 ^^ 750 )}) = 
6420(29) MeV with x^/d.o.f. = 0.92. 

Besides I4, we have also constructed other diquark-antidiquark operators, 
e.g., 


{[f^CXiC]xa[qCc^]a:a + C^]xa[c^Cc]xa] 

to see whether they have good overlap with any resonances as —> rriu- 

However, it turns out that they all behave similar to V 4 . 
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Figure 8 : The mass of the lowest-lying state of M 3 versus the quark mass 
niqa, on the 24^ x 48 lattice at /3 = 6.1. The solid line is the linear ht with 
xVd.o.f. = 0.25. 

5 Summary and Discussions 

In this paper, we have investigated the mass spectra of several interpolating op¬ 
erators (i.e., the hybrid charmonium ccg, the molecular operators Mi, M 2 , and 
M3, and the diquark-antidiquark operator Y 4 ) with = 1 , in quenched 

lattice QCD with exact chiral symmetry. Our results are summarized in Table 
1 , where in each case, the hrst error is statistical, and the second one is our 
estimate of combined systematic uncertainty including those coming from: (i) 
possible plateaus (ht ranges) with /d.o.f. < 1 ; (ii) the uncertainties in the 
strange quark mass and the charm quark mass; (iii) chiral extrapolation (for 
the entries containing u/d quarks); and (iv) hnite size effects (by comparing 
results of two lattice sizes). Note that we cannot estimate the discretization 
error since we have been working with one lattice spacing. Even though lat¬ 
tice QCD with exact chiral symmetry does not have 0{a) and 0{ma) lattice 
artifacts, the 0 {'w?a?) effect might turn out to be not negligible for rrica = 0 . 8 . 

For the hybrid charmonium {eijkC'j^Fjkc), the mass of the lowest-lying state 
only agrees marginally with the mass of Y (4260). Thus it is unlikely to be iden- 
tihed with Y (4260), even though we cannot rule out such a possibility. Never¬ 
theless, we hope to pin down this problem with a more precise determination 
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Figure 9: (a) The time-correlation function C{t) of the lowest-lying state of 
M 3 for niq = rus = 0.08a“^, on the 24^ x 48 lattice at /3 = 6.1. The solid 
line is the hyperbolic-cosine £t for t G [10,16] with /d-o.f. = 1.07. (b) The 
effective mass Meff{t) = \n[C{t)/C{t+ 1)] of C{t) in Fig. 9a. 


1 " ~ ; (SY5YiC){Cy5S)-{Cy5yiS){SY5C) 

24^x48, p = 6.1 
100 configs. 



§ 

TWQCD/l24t48/g100/D5VDA/mq=Q.08 _ 



of the spectrum of the 1 hybrid charmonium in a future publication. 

Evidently, the molecular operator M3 {(q757ic) (cysq) - (cysy^q) (qysc)} 
detects a resonance = 1 ) with mass 4238(31)(57) MeV in the limit 

Viiq —> rriu, which is naturally identified with y(4260). This seems to suggest 
that y(4260) is indeed in the spectrum of QCD, with quark content (cucu). 
Note that we have not studied the excited states of cc, thus we could not rule 
out the possibility that Y (4260) might turn out to be one of the excited states 
of cc, e.g., ip{A^Si) or 7 p{3^Di), even though this is very unlikely in view of 
the widely accepted experimental and theoretical spectrum of cc. 

For the molecular operator Mi rsj {(qyic)(cq) -P (cyiq)(uc)}, and the 
diquark-antidiquark operator I 4 { [q^Cysy.cj [qCysC^j + [q^Cysc] [qCysyf c^j}, 

they also detect states with masses 4350(69) ( 88 ) MeV, and 4267(68) (83) MeV 
respectively, in the limit niq —>■ mu- We suspect that they might be the 
same resonance captured by the molecular operator M3. However, we are 
not sure that these states are resonances since the ratio of spectral weights 
{R = H 20 /H 24 ) for two different lattice volumes with the same lattice spacing 
deviates from one (the criterion for a resonance) with large errors as mq —> mu- 
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Figure 10: The ratio of spectral weights of the lowest-lying state of diquark- 
antidiquark operator I 4 , for 20^ x 40 and 24^ x 48 lattices at /3 = 6.1. The 
upper-horizontal line R = (24/20)^ = 1.728, is the signature of 2-particle 
scattering state, while the lower-horizontal line i? = 1.0 is the signature of a 
resonance. 


It is plausible that such a deviation is due to the quenched artifacts which can 
be evaded if one incorporates internal quark loops. 

Now, in the quenched approximation, our results suggest that Y (4260) has 
a better overlap with the molecular operator M 3 than any other operators 
discussed in this paper. Whether this implies that Y (4260) behaves more 
likely as a HiD molecule than other molecules or diquark-antidiquark meson 
is subjected to further investigations, especially those incorporating dynamical 
quarks. 

Finally, all molecular and diquark-antidiquark operators with quark fields 
(cscs) detect a resonance around 4450±100 MeV, and those with (cccc) detect 
a resonance around 6400 ± 50 MeV. These serve as predictions of lattice QCD 
with exact chiral symmetry. 
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Figure 11: The mass of the lowest-lying state of the diquark-antidiquark op¬ 
erator ^4 versus the quark mass m^a, on the 24^ x 48 lattice at /3 = 6.1. The 
solid line is the linear ht with x^/d.o.f. = 0.49. 
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Figure 12: (a) The time-correlation function C{t) of the lowest-lying state of 
Y 4 for niq = nis = 0.08a“^, on the 24^ x 48 lattice at /3 = 6.1. The solid 
line is the hyperbolic-cosine £t for t G [12,17] with /d.o.f. = 1.03. (b) The 
effective mass Meff{t) = \n[C{t)/C{t+ 1)] of C{t) in Fig. 12a. 
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Figure 13: (a) The time-correlation function C{t) of the lowest-lying state of 
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Operator 

Mass (MeV) 

Resonance/Scattering 


4501(178)(215) 

Resonance 

;^[(u 7 ic)(cu) + (c 7 iu)(uc)] 

4350(69)(88) 

Resonance ? 

^[(s 7 ic)(cs) - 1 - (c 7 is)(sc)] 

4546(30)(61) 

Resonance 

^[(u757ic)(c75u) - (c757iu)(u75c)] 

4238(31)(57) 

Resonance 

;^[(s757ic)(c75s) - (c757is)(s75c)] 

4405(31)(44) 

Resonance 

(c 7 ic)(ss) 

4581(96)(115) 

Resonance 

(c 7 ic)(cc) 

6411(25)(43) 

Resonance 

75 {[u^C'757ic] [ 110750 ^] + [u^075c][u0757fc'^]} 

4267(68)(83) 

Resonance ? 

75 {[ 8 ^ 0757 * 0 ] [ 8 ( 7750 ^] + [s'^C' 75 c][sC' 757 fc^]} 

4449(40)(55) 

Resonance 

75 {(^^<^ 757 *^) ( 0 ( 7750 ^) + (c^( 775 c)(c( 7757 fc^)} 

6420(29)(32) 

Resonance 


Table 1: Mass spectra of hybrid charmonium, molecules, and diquark- 
antidiquark operators with = 1 . 
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